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Abstract 

In a previous paper we found that in the context of the string theory "discretuum" proposed 
by Bousso and Polchinski, the cosmological constant probabihty distribution varies wildly. How- 
ever, the successful anthropic predictions of the cosmological constant depend crucially on the 
assumption of a flat prior distribution. We conjectured that the staggered character of our Bousso- 
Polchinski distribution will arise in any landscape model which generates a dense spectrum of 
low-energy constants from a wide distribution of states in the parameter space of the fundamental 
theory. Here we calculate the volume distribution for A in the simpler Arkani-Hamed-Dimopolous- 
Kachru landscape model, and indeed this conjecture is borne out. 
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I. INTRODUCTION 



While inflationary cosmologists [l|, 0, 0] have long since realized that inflation generically 
gives rise to a multiverse, much more recently string theorists have arrived at a complemen- 
tary world view 0, 0, 0] . Despite a quest to uncover a single unique solution to the Laws 
of Nature, it seems as though string theory admits a vast array of possible solutions. Each 
solution, or vacuum state, represents a possible type of bubble universe, governed by its own 
low-energy laws of physics. 

One can depict each string theory vacuum solution as a local minimum in a multidimen- 
sional potential energy diagram known as the string theory landscape as illustrated in Fig. 
[TJ This landscape of possibilities is expected to have many high-energy metastable false 
vacua which can decay through bubble nucleation0, S, Bubbles of lower-energy vacuum 
can nucleate and expand in the high-energy vacuum background. If the "daughter" vacuum 
has a positive energy density, then inverse transitions are also possible, allowing bubbles of 



high-energy vacuum to nucleate within low-energy vacua [10|, lUj. But if the "daughter" 
vacuum has negative or zero-energy, recycling cannot take place. We will call vacua from 
which new bubbles can nucleate non-terminal, or recyclable vacua, while those which do not 
recycle will be called terminal vacua. This recycling process will populate the multiverse 
with bubbles nested within bubbles of each and every possible type. 

Most of these bubbles will never be hospitable to life. For example, bubbles with large 
positive cosmological constant do not allow for structures such as galaxies or atoms to form 



|12l . |13| . And bubbles with large negative cosmological constant collapse long before life has 
a chance to evolve. However, because the landscape of possibilities is so large, there will also 
be many bubbles which do provide a suitable environment for life to flourish. Obviously we 
live in one of these "friendly" bubbles. 

In the framework of the multiverse, some physical parameters that were once thought of 
as fundamental universal parameters get demoted to local environmental parameters. We 
no longer expect to calculate the "constants" from flrst principles. Instead we are compelled 
to calculate how the "constants" are distributed throughout the multiverse. If we assume we 
are a typical civilization, we should expect to observe values near the peak of the distribution 

3- 

The multiverse paradi gm has led to the successful so-called anthropic prediction of the 



cosmological constant [12, ll3|, [ij, ll5|, ll6|, llTl, ll8|, ll9| . Theoretically we expect the magnitude 



of the cosmological constant^ A ~ 1, but the observed value is Aq ~ 10^^^°. This has been 
one of the biggest problems in theoretical physics^. 

The probability for a randomly picked observer to measure a given value of A can be 
expressed as [13] 

Pobs(A) oc P(A)nobs(A), (1) 

where -P(A) is the volume fraction of regions with a given value of A and nobs (A) is the 
number of observers per unit thermalized volume. 

The factor nobs(A) takes into account selection effects and is sometimes called the an- 
thropic factor. It is in general very difficult to calculate. However, it has been shown 13 



^ Throughout this paper we use reduced Planck units, Mp/Sir = 1, where Mp is the Planck mass. 

^ Furthermore, Aq ^ Pmo where pmo is the present matter density. The smallness of A, and the fact that it 

happens to coincide with the present matter density are collectively known as the cosmological constant 

problems. 
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FIG. 1: The string theory landscape. The vertical axis represents the value of potential energy 
density (or equivalently, the cosmological constant), and the two horizontal directions represent 
two out of hundreds of directions in parameter space, which characterize each vacuum solution. 
Each valley represents a metastable vacuum solution which may have negative, zero or positive 
cosmological constant. 



that the function nobs(^) is only substantially different from zero in a tiny window of width 

AA^ ~ lOOAo ~ 10-^^^ (2) 

around A = 0. AA^ is sometimes called the Weinberg window or the anthropic range. 

The volume factor P(A) depends on the dynamics of eternal inflation and on the under- 
lying fundamental theory. However, it has been argued 2^, 21 1 that it should be accurately 
approximated by a flat distribution. 



P(A) ^ const 



(3) 



because the anthropic range ([2]) is vastly less than the expected Planck scale range of varia- 
tion of A. A smooth function varying on this large characteristic scale will be nearly constant 
within the minute anthropic interval. 

A tiny non-zero value for A was predicted jl, [l2|, 13, 14] when the theoretical vogue was 
to believe that a deep symmetry forced the cosmological constant to be zero. One should 
keep in mind, however, that the successful anthropic prediction for A depends critically on 
the assumption of a fiat volume distribution ([3]). 

In jlll we used the new prescription introduced in [13] to calculate bubble abundances 
in an eternally inflating spacetime to actually calculate the volume distribution for the 
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cosmological constant A in the context of the Bousso-Polchinski (hereafter BP) landscape 
model. We found that the resulting distribution has a staggered appearance, in conflict with 
the heuristically expected flat distribution. 

One might think that the staggered distribution is a feature of the BP model. However, in 
this paper we calculate the volume distribution for A in a simpler landscape model proposed 
by Arkani-Hamed, Dimopolous and Kachru 2J] (hereafter called the ADK model) and once 
again find a wildly varying distribution for A. 

The outline of this paper is as follows: In section [TTl we will describe the ADK landscape 
model. In section HTTl we use the new prescription 23| for calculating probabilities to calculate 
the volume distribution for A in the ADK model. We will conclude with a discussion in 
section [IVl 



II. THE ADK LANDSCAPE 



Following Ref. [2^ we consider a single scalar field with a general quartic potential. 
Also we assume that the theory has two minima at 0± with vacuum energies V±, and we 
take V+ > V-. Thus V+ (VI) represents the energy of the false vacuum (true vacuum) at 
0+ (0_)(see Fig. [2]). We can label the vacua by ?7 = ±1, and define 

n = K. + r/Kiif (4) 

with 

v:„ = 

2 

and 

Kiif=^(n-^-) (6) 

Now consider a theory with J scalar fields (pa, cl = ^,---,J, each having independent 
quartic potentials K(0a) such that the potential is the sum of J independent potentials 



Vav = kv+ + V_) (5) 



V = J2U<Pa) (7) 



a=l 



This theory represents a landscape of 2"^ vacua labeled by {rj} = {rji, ...,rij} with rja = ±1. 
The cosmological constant is given by 

J 

A{r,} =A + Y,VaVma (8) 

a=l 

where 

J 

A = 5^K.a = JK. (9) 

a=l 

If we consider a specific bubble which is completely specified by {rji, rja, ■■■'r]j}, the field 
configuration in an inflating region can change to {rji, rja ± 2, ...//j} when one of the fields 
(pa tunnels to its other minimum through bubble nucleation. Thus, as in the BP model, 
the universe can start off with an arbitrary large cosmological constant, and then diffuse 
through the ADK landscape of possible vacua as bubbles nucleate one within the other. 
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V(cp) 




FIG. 2: Potential with two metastable minima. 

ADK asked the question, does a landscape with 2"^ vacua guarantee that we can solve the 
cosmological constant problem? In the ADK model the histogram of the number of vacua 
per bin of A is well approximated by a Gaussian distribution because A is the sum of many 
independent components. If A = is on the tail of the Gaussian the cosmological constant 
will not scan around 0, but if it is near the peak we can expect to find a dense enough 
spectrum of vacua to account for Aq. Thus ADK concluded that for the landscape to solve 
the cosmological constant problem, either a tiny cosmological constant can accidentally arise 
on the tail of the Gaussian^ or the Gaussian must be peaked close enough to A = 0. 

We wish to point out that even if the Gaussian is peaked around A = 0, yielding a 
sufficiently dense spectrum of A in the anthropic range AA^, this is still not good enough 
to validate the anthropic resolution of the cosmological constant problems. It could be 
that the probabilities of these anthropic vacua differ dramatically (we have learned from 
our calculation in the BP model [12] that the probabilities tend to span many orders of 
magnitude) with one or two dominating the distribution. 

We will now calculate the volume distribution P{A) for the ADK model. 



^ Of course if this were the case then the anthropic explanation for Aq would not be applicable - we simply 
land up at Aq accidentally. 
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III. PROBABILITIES IN THE ADK LANDSCAPE 



In this section we study the volume distribution for A in the ADK landscape using the 
prescription of 23|] which we outline below. We expect to find a wildly varying distribution 
like the one found for the BP model studied in 2^ . 



A. Summary of probability prescription 

Suppose we have a theory with a discrete set of vacua, labeled by index j, and having 
cosmological constants Aj. The volume distribution is given by jl^ 

P.cxp.Zj, (10) 

where pj is the relative abundance of bubbles of type j and Zj is (roughly) the amount of 
slow-roll inflationary expansion inside the bubble after nucleation (so that Zj is the volume 
slow-roll expansion factor). 

The bubble abundances Pi can be related to the comoving volume fractions fj{t) which 
obey the evolution equation ll[ 

i 

where the first term on the right-hand side accounts for loss of comoving volume due to 
bubbles of type i nucleating within those of type j, and the second term reflects the increase 
of comoving volume due to nucleation of type-j bubbles within type-i bubbles. 

The transition rate is defined as the probability per unit time for an observer who is 
currently in vacuum j to find herself in vacuum i and is given by 

Ayr 

^^,=r.,-Hr\ (12) 

where Tij is the bubble nucleation rate per unit physical spacetime volume (same as Xij in 
0) and 

Hj = (A,/3)^/^ (13) 

is the expansion rate in vacuum j. 

Eq. (fTTl) can be written in a vector form, 

§^Mf, (14) 

where f(t) = {fj{t)} and 

Mij = Kij Sij ^ ^ Uri- (1^) 

r 

Let's label the eigenvalue (of the transition matrix M) which has the smallest negative 
real part (in magnitude), —q, and the corresponding eigenvector, s. Then it can be shown 
that the bubble abundances pj are given by 

Pj OC ^ H^KjaSa- (16) 
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where the summation is over all recyclable vacua which can directly tunnel to j. 

The problem of calculating pj has thus been reduced to finding the dominant eigenvalue 
q and the corresponding eigenvector s of the transition matrix M. To calculate the elements 
of M we need to calculate the bubble nucleation rates specific to the landscape model we 
are studying. 



B. Nucleation rates in the ADK landscape 

Transitions between neighboring vacua, which change one of the integers rja by ±2 can 
occur through bubble nucleation. The bubbles are bounded by thin branes, with tension Ta- 
By analogy with the BP model^ we will take the brane tension in the ADK model to be 

Ta = V^AA; = y^W^a (17) 

Transitions with multiple brane nucleation, in which several r/a are changed at once, are 
likely to be strongly suppressed [1^, and we shall disregard them here. 

The bubble nucleation rate Tij per unit spacetime volume can be expressed as 0] 

1^, = Ai,exp-^»^ (18) 

with 

B,, = h,-S, (19) 
Here, lij is the Coleman-DeLuccia instanton action and 

= (20) 

is the background Euclidean action of de Sitter space with expansion rate 

= v/V3- (21) 

In the relevant case of a thin-wall bubble, the instanton action lij has been calculated in 
Refs. 0, 0] • It depends on the values of A inside and outside the bubble and on the brane 
tension r. 

Let us first consider a bubble which changes one rja from rja = +1 to T]a = —1. The 
resulting change in the cosmological constant is given by 

|AA,| = 2\/dif,a (22) 

and the exponent in the tunneling rate (fT8|) can be expressed as 

5ai = 5fr^"^V(x,y). (23) 

j^fiatspace spacc bouucc actiou, 

97 2 4 

Tjflatspace ^ ' 'a fOA^ 

«i " 2 IAAJ3- ^^'^^ 



More specifically for the case of Ua = 1 where Ua is a flux quantum in the BP model. 
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With the aid of Eqs. (fT7|) . (l22i) it can be expressed as 



977^2 

ry flat-space ^ i 'i ("25^ 



The gravitational correction factor r{x, y) is given by 



^ 2[(l + x,)-(l + 2x, + x^).] 
x2(y2_ i)(i + 2x?/ + a;2)5 



with the dimensionless parameters 



-^^ = 1 (27) 



and 



4|AAJ 4 
2A 



I A A. 



(28) 



where A is the background value prior to nucleation. 
The prefactors Aij in f|T8l) can be estimated as 2^ 



A., ~ 1 (29) 

If the vacuum {rji...rja-i, ria — 2, rja+i-...} still has a positive energy density, then an upward 
transition from {r]i...r]a-i,i]a - 2,'r]a+i....} to {r]i...'r]a-i,i]a,i]a+i....} is also possible. The 
corresponding transition rate is characterized by the same instanton action and the same 
prefactor [lo|, and it follows from Eqs. ( fTSl) . ( fT9|) and ( !2TI) that the upward and downward 
nucleation rates are related by 



Fjj — Fjj exp 



-2471^ 



1 1 

Ai Aj 



(30) 



where Aj > Aj. As expected the transition rate from 77^ = — 1 up to r^^ = +1 is suppressed 
relative to that from T]a = +1 down to rja = —1. The closer we are to Aj = 0, the more 
suppressed are the upward transitions i — > j relative to the downward ones. 

We will now investigate the dependence of the tunneling exponent Bai on the parameters 
of the model in the limits of small and large A. For A <^ |AAa|, we have y ~ —1, and Eq. 
gives 

r(y-^-l) = 16. (31) 

The inclusion of gravity increases the tunneling exponent causing a suppression of the nu- 
cleation rate. 

In the opposite limit, A ^ |AAa|, y ^ 1, 

r{y > 1) ^ V2{xy)-^/^, (32) 

and 

Ba^^ 277r'y^J—j . (33) 
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FIG. 3: The spectrum of vacua for a J = 10 ADK landscape with parameters given in Eq. ()35|) 

For large values of A, r -C 1, so the nucleation rate is enhanced. The tunneling action 
must always be large enough to justify the use of the semi-classical approximation: Bai ^ 1, 
or 

A « 30V;;£- (34) 

If Vdif,a and A are changed simultaneously, keeping the ratios A/Vdif,a fixed, then x and y 
do not change, and it is clear from Eqs. ( l23i) .( 1251) .( 126]) that the nucleation exponents scale 
as Bij oc A~^. This shows that bubble nucleation rates are strongly suppressed when the 
energy scales of Vdif,a and A are well below the Planck scale. 

C. Bubble abundances in the ADK model 

We will calculate the bubble abundances for a J = 10 ADK model, containing 2^° vacua, 
with parameter values: 

Vdif,a = {0.0514, 0.0814, 0.0885, 0.1081, 0.1378, 0.1475, 0.1790, 0.2226, 0.2467, 0.2523} (35) 

and A = 0. A histogram of the number of vacua vs. A for this model is given in Fig. [3l 

In direct analogy with our calculation of probabilities for the BP landscape j22|], we 
resort to perturbative techniques, where we use the smallness of upward transitions as a 
small parameter. The slowest vacuum to decay was singled out as the dominant vacuum a^,. 
To zero'th order in perturbation theory (hereafter PT) the only vacua which aquire non-zero 
probabilities are those that are direct offspring from the dominant vacuum a^,. These vacua 
will have large negative cosmological constants. 

The results for the first order bubble abundance factors pj are shown in FigJH 

The dominant site in Fig. H] has "coordinates" (-1, 1,-1, 1, 1,-1, 1,-1,-1, 1) and has a very 
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FIG. 4: Plot of logiQ{l /pj) vs. Aj for the ADK model with parameters given in ()35p and x 
(see (j27p ). The star marks the dominant vacuum a^:. 
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small^ cosmological constant A* = 0.0019. The five squares represent the vacua which can 
be reached from via one upward jump. Their bubble abundances are so low because A,, is 
so small resulting in very suppressed upward jumps (see Eq. (l30l) ). Each site represented by 
a square can then jump down in 5 ways (excluding jumps back to the dominant site itself) 
to the sites depicted as triangles which can in turn jump down to the circles followed by 
crosses. 

The probabilities pj shown here are more suppressed than the BP results in 2^ because 
A* happens to be much smaller. This is simply a consequence of the different parameters. 
Also, unlike the BP results, it appears as though vacua which result from downward jumps 
from a given up jump, have almost the same probability. This "flatness" is fictitious - there 
are actually a few orders of magnitude difference amongst these vacua which is hard to see 
graphically because of the scale. However, overall the "staggered" nature of the distribution 
is evident and the similarity to the BP model is clear. 

In addition to the bubble abundance factor pj, the volume distribution flTUl) includes the 
slow-roll expansion factor Zj. There is no reason to expect the expansion factor to tame 
the wildly varying bubble abundance distribution, and thus we conclude that the volume 
distribution for A will have the same form as that calculated for the bubble abundances. 



IV. CONCLUSIONS AND DISCUSSION 



In |22| we found that for the Bousso-Polchinski string theory landscape, the cosmological 
constant probability distribution varies dramatically for vacua which have close values of A. 



Small relative to the values in the spectrum of our toy model. 
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We have shown that this behavior persists in the case of the Arkani-Hamed-Dimopolous- 
Kachru landscape model. 

This result was expected. Our probability prescription picks a dominant vacuum and all 
other vacua are reached from it via a sequence of upward and downward jumps. To zero'th 
order in PT only the progeny of the dominant vacuum have non-zero probabilities. These 
vacua were shown to have large^ negative A. To the first order in PT only a small subset 
of vacua related to the dominant site via one upward jump and any number of downward 
jumps gain non-zero probabilities. The probabilities of these vacua are proportional to the 
tunneling transition rates of the jumps. The tunneling transition rates have an exponential 
dependence on the parameters of the theory and consequently the probabilities span many 
orders of magnitude, in both landscape models considered. 

Furthermore, it is exceedingly unlikely that one of these first order vacua should be in the 
anthropic range - we just don't have enough of them. So what happens if we go to second 
order in PT? 

Calculations indicate that vacua which can be reached via two upward jumps and sub- 
sequent downward jumps would gain some tiny probabilities. There would be many more 
vacua which can be reached via paths including two upward jumps instead of only one, but 
we would still need to consider higher and higher orders of PT before a sufficient fraction 
of the theory's vacua can be infused with probability. Going to higher orders is technically 
prohibitive. 

Although the distributions we have calculated do not give a fiat distribution to first order, 
we cannot conclude that the anthropic prediction of the cosmological constant was a fluke. 
The distributions we are able to calculate are simply the tip of the iceberg. It is still entirely 
possible that vacua in the anthropic range are smoothly distributed. 

So how do we proceed from here? How do we look beyond the first order results we have 
found for a tiny subset of vacua in our landscape? Currently we do not have a definitive 
answer to this question. But work is underway to try to elucidate what essential features 
of a given landscape model will guarantee that (once we have many vacua in the anthropic 
range) a sufficient number of the most probable anthropic vacua will have close enough 
probabilities to ensure that the distribution can be considered to be smooth. 
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